Variation in the gene coding for the M5 muscarinic receptor (CHRM5) influences cigarette dose but is not associated with dependence to drugs of addiction : evidence from a prospective population based cohort study of young adults

Background: The mesolimbic structures of the brain are important in the anticipation and perception of reward. Moreover, many drugs of addiction elicit their response in these structures. The M5 muscarinic receptor (M5R) is expressed in dopamine-containing neurones of the substantia nigra pars compacta and ventral tegmental area, and regulates the release of mesolimbic dopamine. Mice lacking M5R show a substantial reduction in both reward and withdrawal responses to morphine and cocaine. The CHRM5, the gene that codes for the M5R, is a strong biological candidate for a role in human addiction. We screened the coding and core promoter sequences of CHRM5 using denaturing high performance liquid chromatography to identify common polymorphisms. Additional polymorphisms within the coding and core promoter regions that were identified through dbSNP were validated in the test population. We investigated whether these polymorphisms influence substance dependence and dose in a cohort of 1947 young Australians.

Results: Analysis was performed on 815 participants of European ancestry who were interviewed at wave 8 of the cohort study and provided DNA. We observed a 26.8% increase in cigarette consumption in carriers of the rs7162140 T-allele, equating to 20.1 cigarettes per week (p=0.01). Carriers of the rs7162140 T-allele were also found to have nearly a 3-fold increased risk of developing cannabis dependence (OR=2.9 (95%CI 1.1-7.4); p=0.03).

Conclusion: Our data suggest that variation within the CHRM5 locus may play an important role in tobacco and cannabis but not alcohol addiction in European ancestry populations. This is the first study to show an association between CHRM5 and substance use in humans. These data support the further investigation of this gene as a risk factor in substance use and dependence.

Multiple origins of the endemic Australian Macrobrachium (Descapoda: Palaemonidae) based on 16S rRNA mitochondrial sequences

The evolutionary relationships of the freshwater prawn genus Macrobrachium are obscure. Members of this genus are widely distributed across tropical and subtropical regions. The phylogenetic relationships among the seven endemic and six non-endemic Australian Macrobrachium, along with five non-Australian species, were inferred from the mitochondrial 16S rRNA gene sequences. Methods of analysis yielded phylogenetic trees of differing topologies; however, none supported a monophyletic origin for endemic Australian Macrobrachium. Enforced monophyly of a single origin of endemic Macrobrachium was statistically tested and rejected. These results support the view that the endemic Australian Macrobrachium arose from multiple origins. Previous biogeographical hypotheses related to the radiation of Macrobrachium into Australia are re-examined in the context of these results.

MammoSapiens: eResearch of the lactation program.

Delivering bioinformatics power to life science researchers inevitably runs into problems of limited computing resources in the context of exponentially increasing data sources, access time, costs, lack of skills and, rapidly evolving technology and software tools with poorly defined standards. In this context the development of e-facilities to best enable collaborative research often needs to be customized to specific project applications in close cooperation with the experimentalist users and, to be concerned with the storage and management of results to allow more consistency and traceability of e-results on a broad access data mining platform.

Here we showcase an internet based eResearch platform using the PHP/MySQL paradigm for the collaborative, integrative and comparative analysis of lactation related gene sequences and gene expression experiments to support lactation research. We also illustrate how these resources are used, how they enable research by allowing meta-analysis of data and results and, how the bottom-up development of customized eResearch components can lead to the production of more generic functional software tools and eResearch environments for deployment to a larger number of biological research users working on other bio-systems.

Evolutionary relationships among bluetongue and related orbiviruses

Polymerase chain reaction (PCR) sequencing of specific viral gene segments was used to investigate the phylogenetic relationships among the orbiviruses. Sequence comparisons of the bluetongue virus (BTV) RNA3 from different regions of the world (North America, South Africa, India, Indonesian, Malaysia, Australia and the Caribbean region) showed that geographic separation had resulted in significant divergence, consistent with the evolution of distinct viral populations. There were at least 3 topotypes (Gould, 1987); the Australasian, African - American and another topotype represented by BTV 15 isolated in Australia in 1986. The topotypes of BTV had RNA3 nucleotide sequences that differed by approximately 20 per cent. Analysis of BTV-specific gene segments from animal and insect specimens showed that bluetongue viruses had entered northern Australia from South East Asia, possibly by wind-borne vectors. Nucleotide sequence comparisons were used to show the close genetic relationship between BTV 2 (Ona-A strain) from Florida and BTV 12 from Jamaica, and to investigate the reassortment of BTV genome segments in nature. The mutation rates of the BTV RNA2 and RNA3 segments were estimated to be of the order of 10(-4) nucleotide changes/site/year, similar in magnitude to that reported for other RNA viruses.

Identification and characterization of rhizospheric microbial diversity by 16S ribosomal RNA gene sequencing

In the present study, samples of rhizosphere and root nodules were collected from different areas of Pakistan to isolate plant growth promoting rhizobacteria. Identification of bacterial isolates was made by 16S rRNA gene sequence analysis and taxonomical confirmation on EzTaxon Server. The identified bacterial strains were belonged to 5 genera i.e. Ensifer, Bacillus, Pseudomona, Leclercia and Rhizobium. Phylogenetic analysis inferred from 16S rRNA gene sequences showed the evolutionary relationship of bacterial strains with the respective genera. Based on phylogenetic analysis, some candidate novel species were also identified. The bacterial strains were also characterized for morphological, physiological, biochemical tests and glucose dehydrogenase (gdh) gene that involved in the phosphate solublization using cofactor pyrroloquinolone quinone (PQQ). Seven rhizoshperic and 3 root nodulating stains are positive for gdh gene. Furthermore, this study confirms a novel association between microbes and their hosts like field grown crops, leguminous and non-leguminous plants. It was concluded that a diverse group of bacterial population exist in the rhizosphere and root nodules that might be useful in evaluating the mechanisms behind plant microbial interactions and strains QAU-63 and QAU-68 have sequence similarity of 97 and 95% which might be declared as novel after further taxonomic characterization.

A preliminary study of 16S rRNA sequence variation in Australia Macrobrachium shrimps (Palaemonidae: Decapoda) reveals inconsistencies in their current classification

The systematic relationships among Australian palaemonid shrimps have been the subject of speculation for some time. A preliminary phylogenetic study was undertaken to clarify the relationships of five species, Macrobrachium intermedium (Stimpson), M. australiense (Holthuis), M. atactum (Riek), M. rosenbergii (de Man) and Palaemon serenus (Heller), using 16S rRNA mitochondrial gene sequences. Phylogenetic analyses indicated inconsistencies with the current classification in two respects. First, M. intermedium formed a very well-supported clade with P. serenus distinct from M. australiense, M. atactum and M. rosenbergii. Second, the two species from inland Australia, M. australiense and M. atactum, showed a high level of genetic similarity over a substantial geographic range, suggesting that they may represent conspecific populations. The taxonomic and biogeographic implications of these findings for Macrobrachium in Australia are discussed.

Strategic Australian biotechnology opportunities

Discusses the biotechnology opportunities in Australia as of April 2005. Specific core activities of Australian biotechnology companies; Number of biotechnology companies listed on the Australian Stock Exchange on September 30, 2004; Mutual recognition arrangements that exist with Australia and key international regulatory authorities.

Identification of novel genes differentially expressed in haemopoietic progenitor cells

The biochemical and molecular processes that maintain the stem cell pool, and govern the proliferation and differentiation of haemopoietic stem/progenitor cells (HSPCs) have been widely investigated but are incompletely understood. The purpose of this study was to identify and characterise novel genes that may play a part in regulating the mechanisms that control the proliferation, differentiation and self-renewal of human HSPCs. Reverse transcription differential display polymerase chain reaction (dd-PCR) was used to identify differences in gene expression between a HSPC population defined by expression of the CD34 phenotype, and the more mature CD34 depleted populations. A total of 6 differentially expressed complementary deoxyribonucleic acid (cDNA) sequences were identified. Four of these transcripts were homologous to well characterised genes, while two (band 1 and band 20) were homologous to unknown and uncharacterised partial gene sequences on the GenBank database and were thus chosen for further investigation. The partial cDNA sequences for band 1 and band 20 were designated ORP-3 and MERP-1 (respectively) due to homologies with other well-characterised gene families. Differential expression of the ORP-3 and MERP-1 genes was confirmed using Taqman™ real-time polymerase chain reaction (PCR) with 3 - 4-fold and 4-10 -fold higher levels in the CD34+ fractions of haemopoietic cells compared to CD34- populations respectively. Additionally, expression of both these genes was down regulated with proliferation and differentiation of CD34+ cells further confirming higher expression in a less differentiated subset of haemopoietic cells. The full coding sequences of ORP-3 and MERP-1 were elucidated using bioinformatics, rapid amplification of cDNA ends (RACE) and PCR amplification. The MERP-1 cDNA is 2600 nucleotides (nt) long, and localizes by bioinformatics to chromosome 7.. It consists of three exons and 2 introns spanning an entire length of 31.4 kilobases (kb). The MERP-1 open reading frame (ORF) codes for a putative 344 amino acid (aa) type II transmembrane protein with an extracellular C-terminal ependymin like-domain and an intracellular N-terminal sequence with significant homology to the cytoplasmic domains of members of the protocadherin family of transmembrane glycoproteins. Ependymins and protocadherins are well-characterised calcium-dependant cell adhesion glycoproteins. Although the function of MERP-1 remains to be elucidated, it is possible that MERP-1 like its homologues plays a role in calcium dependent cell adhesion. Differential expression of the MERP-1 gene in haemopoietic cells suggests a role in haemopoietic stem cell proliferation and differentiation, however, its broad tissue distribution implies that it may also play a role in many cell types. Characterization of the MERP-1 protein is required to elucidate these possible roles. The ORP-3 cDNA is 6631nt long, and localizes by bioinformatics to chromosome 7pl5-p21. It consists of 23 exons and 22 introns spanning an entire length of 183.5kb. The ORP-3 ORF codes for a putative 887aa protein which displays the consensus sequence for a highly conserved oxysterol-binding domain. Other well-characterised proteins expressing these domains have been demonstrated to bind oxysterols (OS) in a dose dependant fashion. OS are hydroxylated derivatives of cholesterol Their biological activities include inhibition of cholesterol biosynthesis and cell proliferation in a variety of cell types, including haemopoietic cells. Differential expression of the ORP-3 gene in haemopoietic cells suggests a possible role in the transduction of OS effects on haemopoietic cells, however, its broad tissue distribution implies that it may also play a role in many cell types. Further investigation of ORP-3 gene expression demonstrates a significant correlation with CD34+ sample purity, and 2-fold higher expression in a population of haemopoietic cells defined by the CD34+38- phenotype compared to more mature CD34+38+ cells. This finding, taken together with the previous observation of down-regulation of ORP-3 expression with proliferation and differentiation of CD34+ cells, indicates that ORP-3 expression may be higher in a less differentiated subset of cells with a higher proliferative capacity. This hypothesis is supported by the observation that expression of the ORP-3 gene is approximately 2-fold lower in differentiated HL60 promyelocytic cells compared to control, undifferentiated cells. ORP-3 expression in HL60 cells during normal culture conditions was also found to vary with expression positively correlated with cell number. This indicates a possible cell cycle effect on ORP-3 gene expression with levels highest when cell density, and therefore the percentage of cells in G(0)/G(1) phase of the cell cycle is highest. This observation also correlates with the observation of higher ORP-3 expression in CD34+38-cells, and in CD34+ and HL60 cells undergoing OS induced and camptothecin induced apoptosis that is preceded by cell cycle arrest at G(0)/G(1). Expression of the ORP-3 gene in CD34+ HSPCs from UCB was significantly decreased to approximately half the levels observed in control cells after 24 hours incubation in transforming growth factor beta-1 (TGFâl). As ≥90% of these cells are stimulated into cell cycle entry by TGFâl, this observation further supports the hypothesis that ORP-3 expression is highest when cells reside in the G(0)/G(1) phase of the cell cycle. Data obtained from investigation of ORP-3 gene expression in synchronised HL60 cells however does not support nor disprove this hypothesis. Culture of CD34+ enriched HSPCs and HL60 cells with 25-OHC significantly increased ORP-3 gene expression to approximately 1.5 times control levels. However, as 25-OHC treatment also increased the percentage of apoptotic cells in these experiments, it is not valid to make any conclusions regarding the regulation of ORP-3 gene expression by OS. Indeed, the observation that camptothecin induced apoptosis also increased ORP-3 gene expression in HL60 cells raises the possibility that up-regulation of ORP-3 gene expression is also associated with apoptosis, Taken together, expression of the ORP-3 gene appears to be regulated by differentiation and apoptosis of haemopoietic progenitors, and may also be positively associated with proliferative and G(0)/G(1) cell cycle status indicating a possible role in all of these processes. Given the important regulatory role of apoptosis in haemopoiesis and differential expression of the ORP-3 gene in haemopoietic progenitors, final investigations were conducted to examine the effects OS on human HSPCs. Granulocyte/macrophage colony forming units (CFU-GM) generated from human bone marrow (ABM) and umbilical cord blood (UCB) were grown in the presence of varying concentrations of three different OS - 7keto-cholesterol (7K-C), 7beta-hydroxycholesterol (7p-OHC) and 25-hydroxycholesterol (25-OHC). Similarly, the effect of OS on HL60 and CD34+ cells was investigated using annexin-V staining and flow cytometry to measure apoptosis. Reduction of nitroblue tetrazolium (NBT) was used to assess differentiative status of HL60 cells. CFU-GM from ABM and HL60 growth was inhibited by all three OS tested, with 25-OHC being the most potent. 25-OHC inhibited ≥50% of bone marrow CFU-GM and ≥95% of HL60 cell growth at a level of 1 ug/ml. Compared to UCB, CFU-GM derived from ABM were more sensitive to the effects of all OS tested. Only 25-OHC and 7(5-OHC significantly inhibited growth of UCB derived CFU-GM. OS treatment increased the number of annexin-V CD34+ cells and NBT positive HL60 cells indicating that OS inhibition of CFU-GM and HL60 cell growth can be attributed to induction of apoptosis and differentiation. From these studies, it can be concluded that dd-PCR is an excellent tool for the discovery of novel genes expressed in human HSPCs. Characterisation of the proteins encoded by the novel genes ORP-3 and MERP-1 may reveal a regulatory role for these genes in haemopoiesis. Finally, investigations into the effects of OS on haemopoietic progenitor cells has revealed that OS are a new class of inhibitors of HSPC proliferation of potential relevance in vivo and in vitro.

Accurately measuring recombination between closely related HIV-1 genomes

Retroviral recombination is thought to play an important role in the generation of immune escape and multiple drug resistance by shuffling pre-existing mutations in the viral population. Current estimates of HIV-1 recombination rates are derived from measurements within reporter gene sequences or genetically divergent HIV sequences. These measurements do not mimic the recombination occurring in vivo, between closely related genomes. Additionally, the methods used to measure recombination make a variety of assumptions about the underlying process, and often fail to account adequately for issues such as co-infection of cells or the possibility of multiple template switches between recombination sites. We have developed a HIV-1 marker system by making a small number of codon modifications in gag which allow recombination to be measured over various lengths between closely related viral genomes. We have developed statistical tools to measure recombination rates that can compensate for the possibility of multiple template switches. Our results show that when multiple template switches are ignored the error is substantial, particularly when recombination rates are high, or the genomic distance is large. We demonstrate that this system is applicable to other studies to accurately measure the recombination rate and show that recombination does not occur randomly within the HIV genome.

Interspecific variations in the gastrointestinal microbiota in penguins

Despite the enormous amount of data available on the importance of the gastrointestinal (GI) microbiota in vertebrate (especially mammals), information on the GI microbiota of seabirds remains incomplete. As with many seabirds, penguins have a unique digestive physiology that enables them to store large reserves of adipose tissue, protein, and lipids. This study used quantitative real-time polymerase chain reaction (qPCR) and 16S rRNA gene pyrosequencing to characterize the interspecific variations of the GI microbiota of four penguin species: the king, gentoo, macaroni, and little penguin. The qPCR results indicated that there were significant differences in the abundance of the major phyla Firmicutes, Bacteroides, Actinobacteria, and Proteobacteria. A total of 132,340, 18,336, 6324, and 4826 near full-length 16S rRNA gene sequences were amplified from fecal samples collected from king, gentoo, macaroni, and little penguins, respectively. A total of 13 phyla were identified with Firmicutes, Bacteroidetes, Proteobacteria, and Fusobacteria dominating the composition; however, there were major differences in the relative abundance of the phyla. In addition, this study documented the presence of known human pathogens, such as Campylobacter, Helicobacter, Prevotella, Veillonella, Erysipelotrichaceae, Neisseria, and Mycoplasma. However, their role in disease in penguins remains unknown. To our knowledge, this is the first study to provide an in-depth investigation of the GI microbiota of penguins.

Phylogenetic relationships among nine scallop species (Bivalvia : Pectinidae) inferred from nucleotide sequences of one mitochondrial and three nuclear gene regions

Current knowledge of the evolutionary relationships among scallop species (Mollusca: Bivalvia: Pectinidae) in the Indo-Pacific region is rather scanty. To enhance the understanding of the relationships within this group, phylogenies of nine species of scallops with the majority from coastal regions of Thailand, were reconstructed by maximum parsimony, maximum likelihood, and Bayesian methods using sequences of the 16S rRNA of the mitochondrial genome, and a fragment containing the ITS1, 5.8S and ITS2 genes of the nuclear DNA. The trees that resulted from the three methods of analysis were topologically identical, however, gained different levels of support at some nodes. Nine species were clustered into two major clades, corresponding to two subfamilies (Pectininae and Chlamydinae) of the three currently recognized subfamilies within Pectinidae. Overall, the relationships reported herein are mostly in accordance with the previous molecular studies that used sequences of the mtDNA cytochrome oxidase subunit I, and the classification system based on microsculpture of shell features and morphological characteristics of juveniles. Levels of divergences were different among genes (i.e., the 5.8S gene showed the lowest levels of nucleotide divergence at all levels, whereas the 16S rRNA showed the highest level of variation within species, and ITS2 gene revealed the highest level of divergence at higher levels).

De novo assembly and characterization of the invasive Northern Pacific Seastar transcriptome

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI's non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

Utility of mitochondrial DNA sequences from four gene regions for systematic studies of Australian freshwater crayfish of the Genus Cherax (Decapoda: Parastacidae)

One of the most important requirements for systematic and phylogenetic studies is the identification of gene regions with the appropriate level of variation for the question of interest. Molecular phylogenetic and systematic studies of freshwater crayfish have made use of DNA sequences mainly from ribosomal genes, especially the 16S rRNA gene region. Thus, little information is available on other potentially useful mitochondrial gene regions for systematic studies in these animals. In this study, we look at nucleotide variation and phylogenetic relations within and between four species of freshwater crayfish of the genus Cherax from the southwest of Western Australia using four fragments amplified from the 16S rRNA, 12S rRNA, Cytochrome Oxidase I (COI), and Cytochrome b (Cyt b) gene regions. Samples of Engaeus strictifrons, Euastacus bispinosus, and Geocharax falcata were also sequenced for comparative purposes. The size of the fragments varied from 358 bp to 600 bp. Across all samples, the four fragments showed significant phylogenetic signal and showed similar proportions of variable sites (28.81–37.33%). Average divergence within species for the mitochondrial gene regions varied from 1.18% to 4.91%, with the 16S rRNA being the least variable and Cyt b the most variable. Average divergence between species ranged 7.63–15.53%, with 16S rRNA being the least variable and COI the most variable. At the generic level, average divergence ranged 17.21–23.82%. Phylogenetic analyses of the 16S rRNA, 12S rRNA, and COI regions generated four clades consistent with the presence of four species previously identified on the basis of allozyme and morphological studies. The relationships among samples were largely congruent across the data set, although some relationships remained unresolved. Not all samples could be amplified using the Cyt b primers, and some of those that were showed quite anomalous relationships, suggesting that one or more Cyt b pseudogenes were being amplified.

The systematics of freshwater crayfish of the genus Cherax Erichson (Decapoda : Parastacidae) in eastern Australia re-examined using nucleotide sequences from 12S rRNA and 16S rRNA genes

Nucleotide sequence data were used to re-examine systematic relationships and species boundaries within the genus Cherax from eastern Australia. Partial sequences were amplified from the 12S (~365 bp) and 16S (~545 bp) rRNA mitochondrial gene regions. Levels of intra- and inter-specific divergence for Cherax species were very similar between the two gene regions and similar to that reported for other freshwater crayfish for 16S rRNA. Phylogenetic analyses using the combined data provided strong support for a monophyletic group containing 11 eastern Australian species and comprising three well-defined species-groups: the 'C. destructor' group containing three species, the 'C. cairnsensis' group containing four species and the 'C. cuspidatus' group containing two species. Cherax dispar and C. robustus are distinct from all other species and each other. In addition, two northern Australian and a New Guinean species were placed in the 'Astaconephrops' group, which is the sister-group to the eastern Australian Cherax lineage. Several relationships were clarified, including: the status of northern and southern C. cuspidatus as separate species; a close relationship between C. cairnsensis and C. depressus; the validity of C. rotundus and C. setosus as separate species and their close affinities with C. destructor; and the distinctiveness of the northern forms of Cherax. The analysis of the 12S rRNA and 16S rRNA data is highly concordant with the results of previous allozyme studies.